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Abstract 


We generalize the linear cr model in order to develop a chiral-invariant model 
of nuclear structure. The model is natural, and contains not only the usual 
a meson which is the chiral partner of the pion but also a new chiral-singlet 
that is responsible for the medium-range nucleon-nucleon attraction. This 
approach provides significant advantages in terms of its description of nuclear 
matter and finite nuclei in comparison with conventional models based on the 
linear a model. 
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Phenomenological relativistic field theories can provide a satisfactory description of the 
properties of nuclear matter and finite nuclei [1,2]. The fundamental importance of chiral 
symmetry means that it must be incorporated into models of this kind. As a result, there is a 
long history of attempts to generalize the linear a model (LaM) to build nuclear models with 
chiral symmetry. However, it is well known that the LaM, supplemented by a repulsive force 
associated with cu-exchange, leads to bifurcations in the equation of state (EOS) [3,4] and 
that the model usually gives a very large incompressibility, K, and, hence, a very stiff EOS. 
The bifurcation problem may be settled (1) by introducing a a-u coupling in the mean-field 
approximation (MFA) [4] and (2) in the one-loop approximation due to a modified effective 
potential [5]. 

It is, however, difficult to build relativistic mean field phenomenology with manifest 
chiral symmetry based upon the conventional LaM, which contains a scalar meson (namely 
the (7 meson) playing a dual role as the chiral partner of the pion and the mediator of 
the medium-range nucleon-nucleon (NN) attraction [6]. Chiral perturbation theory actually 
forbids the identification of this NN attraction with the exchange of the a. In the LaM, 
two critical constraints are imposed: one is the “Mexican hat” potential, which gives a very 
strong nonlinearity of the scalar field, and the other is the equality between the scalar and 
pion couplings to the nucleon. The latter eventually requires a large value (> 600MeV) 
of the a mass and leads to unrealistic, oscillating nuclear charge densities, level ordering 
and shell closures [6]. Furthermore, as is well known, the predictions of the LaM generally 
involve cancellations among several graphs. One example is the ttN scattering amplitude 
where the a exchange combining with the Born term satisfies the soft pion results, whose 
magnitude is too small. 

An alternative is to adopt a nonlinear realization of chiral symmetry. To avoid the 
unnatural cancellations of the LciM, the a field is eliminated by imposing the constraint 
that the fields he on the chiral circle, (with the pion decay constant), in the 

nonlinear model (NLM). This is all fine for chiral symmetry. However, at mean-field level, 
it is somewhat frustrating for nuclear physics because it is not easy to introduce a scalar- 
isoscalar channel which is responsible for the mid-range NN attraction. When studying only 
elementary processes, the cancellations are just a matter of taking care but when studying 
the nuclear many-body problem in MFA, it is certainly much simpler to introduce a scalar- 
isoscalar meson explicitly. 

In the usual NLM, the chiral radius is fixed to by mere convenience. However, nothing 
prevents us from keeping the fluctuation of the chiral circle as a degree of freedom, and 
identifying it with the meson which yields the mid-range attraction. Chanfray et ah [7] have 
proposed a chiral nuclear model in which the LaM is reformulated in the standard nonlinear 
form as far as the pion field is concerned but where the scalar degree of freedom (called 
6), corresponding to fluctuations along the chiral radius, is maintained. Thus, the model 
explicitly involves the scalar meson for the NN attraction but does not have the unnatural 
cancellations. However, in their calculation the mass of the 6 meson was predicted to be in 
the range 0.8 — 1.0 GeV, which may again be too large to obtain good fits to the properties 
of finite nuclei. 

Combining nonlinear chiral symmetry with the broken scale invariance of QCD is another 
approach to construct a nuclear model [8]. In such a model, the low-energy theorems involv¬ 
ing the QCD trace anomaly of the energy-momentum tensor are assumed to be saturated by 
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a scalar glueball (gluonium) field with a large (above 1 GeV) mass and a light scalar field. 
Thus, two different scales are required in the scalar channel. After integrating out the heavy 
gluonium held, one can construct an effective model which contains the chiral-singlet scalar 
held with a mass of about 400-500 MeV. Since the light scalar meson can be responsible 
for the mid-range NN attraction, the model provides good hts to the bulk properties of h- 
nite nuclei and single-particle spectra. However, this does not involve the scalar held which 
corresponds to the huctuation of the chiral radius. 

Even in the LaM, if a new, light scalar held which phenomenologically simulates the 
mid-range attractive NN force is introduced not as the chiral partner of the pion but as a 
chiral singlet, one may be able to construct a new chiral-invariant model. Here we construct 
such a model and study the role of scalar mesons in symmetric nuclear matter as well as 
hnite nuclei. As discussed above, it may be necessary to include at least two kinds of scalar 
mesons: one the usual a meson with a relatively large mass and the other a new, chiral- 
singlet scalar meson with a light mass. The a meson is the chiral partner of the pion and 
they form an isoscalar-isovector scalar quartet, (n, tt). The huctuations around the stable 
point on the chiral circle are described by them. When the a mass is large, the a held itself 
is relatively reduced in matter and it is hence expected that the “Mexican hat” potential 
producing the undesirable, strong nonlinearity of the a held would also be suppressed. In 
contrast, because of its small mass, the new scalar held is enhanced and it may provide the 
main part of the NN attraction. Our underlying philosophy is that the NN attraction in 
matter is simply dominated by scalar-isoscalar, correlated two-pion exchange and that it is 
this that is being simulated by the light scalar meson in MFA. 

Let us hrst start from the Lagrangian density of the LaM with an explicit symmetry 
breaking term 


^ — -^0 + ^SB, ( 1 ) 

where 

Cq = - go{a + if ■ 7^'y5)]^|J + ^{d^ad'^a + d^f ■ (9^if) - -h if^ - (2) 

and CsB = ca with a and f the a and tt helds, respectively. Next we introdnce a chiral- 
invariant Lagrangian density for the chiral-singlet scalar meson s: 

Cs = + if ■ + ]-{d^,sd^s -mls^), (3) 

jn ^ 

with Qs the strength of the s-N coupling and rUg the s-meson mass. The first term describes 
the interaction between the s meson and nucleon, which is similar to that adopted in Eq.(20) 
of Ref. [9]. In general, F{s) is allowed to be an arbitrary function of the scalar field. 
However, in the present model we choose the simplest one, namely, F{s) = s. Note that, 
as demonstrated in the QMC model [10,11], this complexity is associated with the fact that 
the nucleon has the internal strncture which responds to the mean scalar helds. In principle 
one could also multiply this interaction by an arbitrary function of However, since 

the a held is relatively suppressed because of its large mass, we ignore snch possibilities in 
this initial investigation. We simply add this Lagrangian density with F{s) = s to the LuM. 
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It is also vital to add the repulsive NN interaction due to the cn-meson exchange. It can 
be introduced as a gauge-like boson [4], that is, we add the kinetic energy term for the u 
meson and replace the derivative by a covariant form: = d^ + with the 

cn-N coupling constant. 

As usual, we shift the a held (n —> — a) and eliminate A, v and c in favor of the free 
masses of the nucleon, a and n mesons (M, nia and uIt,, respectively). (Notice that hereafter 
a denotes the shifted, positive mean value.) The free masses are generated by spontaneous 
symmetry breaking. Taking MFA, the total Lagrangian density then reads 

C = - M*{a, s) - + ^{di^ad^a + df^sd^s) - -h m^s^) 

- {a, s)uj^ -V{a), (4) 

where the pion held vanishes, because the nuclear ground state has good parity. The time 
component of the u held is simply denoted by u and the “Mexican hat” potential, V, is 

= (5) 

with R = (mjr/mo-)^ and m^(= 138 MeV) the pion mass. In the chiral limit, R becomes 0. 
The ehective nucleon and u meson masses are respectively given by 


M*{a, s) 


(M - g,s) 



M -goa- gsS + —{goa){gss), 


*2 / \ 2 
(a, s) = 



( 6 ) 

(7) 


Here the free masses and the pion decay constant in vacuum are taken to be M = 939 
MeV, = 783 MeV and f-„ = 93 MeV. Because the pion decay constant is hxed, the 
coupling constants go and g^ are automatically determined through the relations M = f^^go 
and = g^f^. 

This gives the total energy per nucleon (for symmetric nuclear matter) at baryon density 
Pb- 


E 

A 


—— / dk V + k^ + -—( 
{2t[Ypb J 2pB 


2 2 I 2 2\ I 

m^a + m^s ) + 


alpB ^ vjcT) 
Pb 


2mY 


( 8 ) 


There are three parameters to be determined: the scalar meson masses, m„ and and 
the coupling constant Ps- First, let us try to hx ttIo- to be a relatively large mass around 1 
GeV. The remaining parameters, rris and Ps-, are then chosen so as to produce the saturation 
property of symmetric nuclear matter: E/A — M = —15.7 MeV at po = 0.17 fm“^ (po the 
normal nuclear matter density). 

This Lagrangian, however, produces a very small incompressibility, around 100 MeV. As 
seen in Eq.(7), the coupling of the a meson to the u reduces the u mass in matter, while the 
light scalar meson enhances it. Since the u mass is eventually not reduced much at large pb, 
the repulsive force is insufficient and the model cannot produce the correct incompressibility 
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{K = 210 ± 30 MeV). If the effective uj mass is reduced significantly, the EOS becomes stiff 
[12], To cure this defect, let us add a (chiral-invariant) term 


r = 


2 /, 




to the Lagrangian density. This term changes the effective on meson mass as 


m 


*2 

LO 








( 9 ) 


( 10 ) 


The new parameter ( can control m* around po and provide a more realistic incompressibil¬ 
ity. Eq.(lO) also implies that, as in the nucleon, the u meson has internal structure which 
depends on the scalar helds [11]. 

Now we are in a position to show our results for nuclear matter. Varying the a meson 
mass, we determine the three parameters, m*, p* and (, so as to reproduce the saturation 
condition and the correct K. In Table I, we show the nuclear matter properties at po. We 
also present the scalar and vector mean helds as functions of nuclear matter density in Fig.l. 
In our calculation, the a mass varies around 1 GeV. As discussed above, the mass of the 
light scalar meson should be about 400 — 500 MeV to get good hts to properties of hnite 
nuclei. From the table, in the case of m^j = 1.0 — 1.05 GeV, the s meson mass is in the 
desired range. This is not a trivial fact. In contrast, in the case where > 1.1 GeV 
(rUo- < 0.95 GeV), m* seems too large (small). Furthermore, as expected, with increasing 
a mass, the contribution of the a held to the total energy is suppressed (see the table and 
hgure) and the “Mexican hat” potential is correspondingly much reduced. Note that, in the 
chiral model of Boguta [4], we hnd g^a ~ 200 MeV and V/pb — —21 MeV at po. 

The nucleon mass is reduced by about 25% at po. In Eq.(6), the linear terms of the 
scalar helds reduce the mass while the quadratic term increases it. It should be emphasized 
here that the reduction of the nucleon mass depends on the s held as well as the a, where 
the a may be directly related to the change of quark condensate in matter (see discussions 
below), while the s-meson contribution is related to the mid-range NN attractive force. In 
this model, the two diherent origins of the mass reduction are simultaneously present. 

The contribution of the quadratic term in M* has already been studied in detail in the 
QMG model from the point of view of the quark substructure of the nucleon [10,11]. In the 
QMG model, the ehective nucleon mass is approximated as M* k, M — g^s + ^{gssYi where 
the coefficient of the quadratic term, a/2, is known as the scalar polarizability [10,13] (a is 
estimated to be 10“^ MeV“^ [H])- Ihe present case, since numerically g^a ~ 0.5 x ggS for 
the proper a mass (see the table I), the scalar polarizability is of order 1/2M ~ 0.5 x 10“^ 
MeV“^, which is just the same as that in the QMG model. 

We now wish to relate the change of the nucleon mass to the evolution of the quark 
condensate {gqips)) in matter. The quark condensate and its evolution at hnite density 
can be obtained by identifying the symmetry breaking pieces of QGD with that of our 
Lagrangian: —2mqqq = ca = f-^rn^a. This shows that the condensate evolution is driven by 
the mean value of a, the chiral partner of the pion. From this relation and the Gell-Mann- 
Oakes-Renner relation, we get the relative modihcation of the condensate (to the vacuum 
value) at hnite density 
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( 11 ) 


(w(o)) W 

Since both the a and s fields at low density are exclusively governed by the scalar density 
of the nucleon field 

a ~ and s ~ (12) 

mj. 

the s field at low density can be related to the a field 



Using Eqs.(6), (11) and (13), we can numerically relate the nucleon mass at low density to 
the condensate evolution 

Af ^ (ot(pb)) 

M (ijgiO)) 

This rather complicated structure comes from the fact that the mass reduction depends two 
different sources, both the quark condensate and the NN attractive force. 

Next let us show some properties of finite nuclei. As an illustrative example we consider 
the case of ^°Ca. Firstly, instead of po = 0.17 fm“^, we adopt the value of po = 0.155 
fm“^, which is closer to the interior density of ^°®Pb, and again search the parameter set 
for nuclear matter. We then find: = 0.92 GeV, = 404.0 MeV, = 18.07 and 

( = 0.7, which reproduces the saturation condition with K = 218 MeV. We next solve a 
set of coupled nonlinear differential equations for a finite nucleus, that is derived from the 
total Lagrangian density (see Eqs.(4) and (9)). It may be solved by a standard iteration 
procedure [10,11]. This leads to a value for the binding energy per nucleon of Eb/A = —8.28 
MeV (the observed value is —8.45 MeV). 

In Table II, the calculated spectra are summarized. Because of the relatively smaller 
scalar and vector fields in the present model than in Quantum Hadrodynamics (QHD) [1], 
the spin-orbit splittings are smaller. The good agreement in the binding energy per nucleon 
comes at the expense of a reduction in the spin-orbit force. Note that there is a strong 
correlation between the effective nucleon mass and the spin-orbit force. It is expected that 
the inclusion of the exchange contributions (Fock terms) will increase the absolute values 
of the single-particle, scalar and vector potentials [14]. This increase, together with the 
contribution to the spin-orbit force from the Fock term itself, may be expected to improve 
the spin-orbit splitting in finite nuclei. We leave this for a future study. 

The charge density distribution is illustrated in Fig.2. Having solved the coupled dif¬ 
ferential equations, we obtain the point-proton and neutron densities in a nucleus. It is 
then necessary to consider the effect of the nucleon form factor that gives a considerable 
correction to the density distribution. We calculate the charge density by a convolution of 
the point-proton and neutron densities with the proton and neutron charge distributions 
[10,15]. We then obtain the root-mean-square (rms) charge radius {tch = 3.44 fm) and the 
difference between nuclear rms radii for neutrons and protons (r„ ~ = —0.157 fm) - the 

observed values are, respectively, 3.48 fm and 0.05 ± 0.05 fm. 


gsTTia 

,goms 


{QQjpB)) ' 
(9g(0)) , 


(14) 
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As seen in the figure, the calculated charge density distribution is close to the experi¬ 
mental area. We note that there are no strong oscillations, even in the charge density of 
(In conventional nuclear models based on the LaM, unrealistic oscillations in the 
charge density have been reported [6].) In Fig.3, we show the neutron density distribution. 
Again we find reasonable agreement with the experimental data. The results for the charge 
and neutron densities in our model are very similar to those of the other chiral models [8]. 

In general, an effective field theory at low energy contains an infinite number of inter¬ 
action terms, which incorporate the compositeness of hadrons [11]. It is then expected to 
involve numerous couplings which may be nonrenormalizable. To make sensible calculations, 
Manohar and Georgi [18] have proposed a systematic way to manage such complicated, effec¬ 
tive field theories called “naive dimensional analysis” (NDA). NDA gives rules for assigning 
a coefficient of the appropriate size to any interaction term in an effective Lagrangian. Af¬ 
ter extracting the dimensional factors and some appropriate counting factors using NDA, 
the remaining dimensionless coefficients are all assumed to be of order unity. This is the 
so-called naturalness assumption. If naturalness is valid, the effective Lagrangian can be 
truncated at a given order with a reasonable bound on the truncation error for physical 
observables. Then we can control the effective Lagrangian, at least at the tree level. 

Here we use NDA to see whether this model gives natural coefficients. When a generic, 
interaction Lagrangian density is written as [19] 


C 


int. 


ce, 


.mnp 


V>r(r/2)^ Y 


m\n\p\ \ /^A J 



( 15 ) 


the overall coupling constant cemnp is dimensionless and of 0{1) if naturalness holds. Here 
T and r, respectively, stand for a combination of Dirac matrices and isospin operators, and 
A(~ 1 GeV) is a large mass scale for the strong interaction. 

The present model has 10 interaction terms and the result is shown in Table HI. Only 
three coefficients, C 0023 , C 0300 and co 4 oo, are close to 3: C 0023 comes from Eq.(9) while the 
latter two are for the ’’Mexican hat” potential. They all are, however, of order unity and 
almost natural. This model is thus natural as an effective field theory for nuclei. 

In summary, we have developed a chiral nuclear model which is natural, based on the 
LaM. To describe the properties of nuclear matter and finite nuclei at mean-field level, it 
seems to be necessary to introduce two different scales associated with the scalar mesons: 
one is the usual a meson that is the chiral partner of the pion, and the other is a new one 
that is treated as a chiral-singlet and is responsible for the mid-range NN attraction. The 
present model has three parameters, which are determined so as to produce the saturation 
condition at normal nuclear matter density and the proper incompressibility. Then, the 
model automatically predicts two scalar scales, that is, one concerns the chiral symmetry 
(ttIo- ~ 1 GeV) and the other relates to the NN attraction (m^ ~ 400 — 500 MeV). Because 
the a mass is heavy, the contribution of the a field to the total energy is relatively suppressed 
and the “Mexican hat” potential is very much reduced as well. This fact is vital to obtain 
better fits to the properties of finite nuclei. We have demonstrated the single particle spectra, 
charge and neutron densities of ^‘^Ga as well as the properties of symmetric nuclear matter. 
The present approach improves conventional nuclear models based on the LaM. As a next 
step, it would be necessary to extend this approach to include the p meson in order to study 
unstable nuclei and/or dense matter like neutron stars. 
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It would be of great interest to study, from the point of view of quark substructure of 
hadrons, the relation between the dependence of hadron masses on the scalar helds in this 
chiral model and that in the QMC model [10,11]. In particular, if the effective lo mass is 
sufficiently reduced around normal nuclear matter density, it would be possible to form an cu- 
nucleus bound state [20] and such exotic states may provide signihcant information on chiral 
symmetry restoration (see Eqs.(lO) and (11)). Finally, we observe that the a mass predicted 
in this model is very close to the scalar glueball mass assumed in Ref. [8] (although it was 
integrated out and does not appear explicitly). Therefore, it would also be very intriguing 
to study whether there is a deeper connection between the a and such scalar glueball. 

Acknowledgment: This work was supported by DOE contract DE-AC05-84ER40150, under 
which SURA operates Jefferson Laboratory. 
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TABLES 


TABLE I. Properties of symmetric nuclear matter at normal nuclear matter density (po = 0.17 
fm“^). The a mass varies in the range of 0.9 and 1.15 GeV. Note that, in the present model, 
gQ = 101.9 and = 70.89, which together with the pion decay constant (/jr = 93 MeV) give the free 
nucleon and uj meson masses. We take ttItt = 138 MeV. All masses, energies and incompressibility 
are quoted in MeV. 


rria 

nis 

al 

c 

goer 

9sS 

y{er)/pQ 

M*/M 

K 


346.5 

12.28 

0.7 

130 

130 

-6.69 

0.743 

214 


388.0 

17.29 

0.75 

111 

143 

-4.72 

0.747 

216 


441.7 

24.48 

0.8 

96.4 

155 

-3.42 

0.750 

214 


510.7 

35.00 

0.85 

84.4 

164 

-2.55 

0.751 

214 


612.1 

53.08 

0.85 

74.6 

172 

-1.93 

0.752 

212 

1150 

741.9 

81.35 

0.85 

66.5 

179 

-1.50 

0.752 

218 


TABLE II. Single particle energies for ^*^Ca. All energies are quoted in MeV. 


l'Si/2 

lP3/2 

lPl/2 l'^5/2 

l'^3/2 

2'Si/2 

proton 

-35.8 

-24.8 

-23.2 -12.7 

-10.0 

-7.7 

neutron 

-44.8 

-33.2 

-31.7 -20.5 

-17.9 

-15.8 


TABLE III. Interaction terms and corresponding (dimensionless] 
rUfj = 1050 MeV (see Tablel). 

1 coupling constants for 

term 

^£mnp 

value 

ipaip 

CllOO 

0.94 

tpsip 

ClOOl 

0.55 

'4)as'4) 

CllOl 

0.55 


ClOlO 

0.78 

OijP' 

C0120 

1.2 


C0220 

1.2 


C0022 

1.2 


C0023 

3.1 


C0300 

3.3 


C0400 

3.3 
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FIGURES 



Pb ^ Po 

FIG. 1. Scalar and vector mean fields. The dotted and solid curves are, respectively, for the a 
(gocr) and s {gss) fields, while the dot-dashed one is for the uj (gu^uj) field. The a mass is taken to 
be 1050 MeV (see Table I). 



r(fm) 

FIG. 2. Charge density distribution for “^^Ca. The experimental data are denoted by the dashed 
area [16]. 
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FIG. 3. Point-neutron density distribution in The dashed area represents the empirical 

fit to proton scattering data performed in Ref. [17]. 
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